% Scripts to plot EMAG2 grid filled with MF6
% Nov 2013
emag = load('/Users/manojnair/Downloads/emagwithMF61NoFlag.xyz');
load coast
imagesc([-180:180],[-90:90],reshape(emag(:,3),[360*30, 180*30+1])');

caxis([-200 200]);

hold
plot(long,lat,'k')
set(gca,'FontSize', 16);
set(gca,'YDir','normal');
set(gca,'xtick',[],'ytick',[])
colorbar('h','FontSize', 16)
